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Abstract 



We introduce the method of dynamical renormalization group to study re- 
laxation and damping out of equilibrium directly in real time and applied 
it to the study of infrared divergences in scalar QED. This method allows 
a consistent resummation of infrared effects associated with the exchange of 
quasistatic transverse photons and leads to anomalous logarithmic relaxation 
of the form e~ aTt ln [*Ao] for hard momentum charged excitations. This is in 
contrast with the usual quasiparticle interpretation of charged collective exci- 
tations at finite temperature in the sense of exponential relaxation of a narrow 
width resonance for which the width is the imaginary part of the self-energy 
on-shell. In the case of narrow resonances away from thresholds, this approach 
leads to the usual exponential relaxation. The hard thermal loop resumma- 
tion program is incorporated consistently into the dynamical renormalization 
group yielding a picture of relaxation and damping phenomena in a plasma in 
real time that trascends the conceptual limitations of the quasiparticle picture 
and other type of resummation schemes. 

12.38.Mh,11.15.-q;11.15.Bt 
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I. INTRODUCTION 



The possibility of studying experimentally the formation and evolution of the quark 
gluon plasma at RHIC and LHC motivates a deeper understanding of collective excitations 
in ultrarelativistic plasmas (for reviews see |10|)- An important concept in the assessment 
of whether the quark gluon plasma achieves local thermodynamic equilibrium is that of the 
relaxation time scale or damping rate which determines the lifetime of excitations in the 
plasma P-l2|. The concept and definition of the damping rate of a collective excitation 



is associated with a quasiparticle description of these excitations in the plasma and imply 
exponential relaxation. The validity of the quasiparticle concept requires that the lifetime 
must be large compared to the oscillation period of the particular excitation mode. In 
this quasiparticle picture the collective excitations are described as narrow resonances, their 
spectral function is of the Breit-Wigner form and the damping rate is obtained from the width 
of this resonance. For weakly interacting quasiparticles the narrow resonance (quasiparticle) 
approximation is expected to be reliable and the damping rate or lifetime is obtained from 
the imaginary part of the self-energy on the mass shell of the collective excitation [p|-l2|. 

Early attempts to calculate the damping rates of quasiparticles in lowest order pertur- 
bation theory obtained gauge dependent and unphysical results ||13|| . Braaten and Pisarski 
[njj- |TB[ introduced a resummation scheme (the resummation of the hard thermal loops or 
HTL) that incorporates the screening corrections in a gauge invariant manner and render 
finite transport cross sections |T3. These hard thermal loop screening corrections are suffi- 



cient to render finite the damping rate of excitations at rest in the plasma. However these 
screening corrections are not sufficient to cure the infrared divergences in the damping rate 
of charged excitations at non-zero and large (hard) spatial momentum |[20| . The infrared 
divergences arise from the emission and absorption of long-wavelength magnetic (transverse) 
photons or gluons which are not screened by the hard thermal loop corrections [^U] . Whereas 
longitudinal photons (instantaneous Coulomb interaction) are screened at finite temperature 
with a Debye screening mass ~ eT, magnetic photons (transverse) are dynamically screened 
for non zero frequency as a result of Landau damping [13|]- |16| , p!8|j20|] . However, quasistatic 
long-wavelength magnetic photons are not screened in the Abelian theory, and their emis- 
sion and absorption by a fast moving charged particle results in infrared divergences in the 
imaginary part of the self-energy on shell. 

These infrared divergences for charged particles are not specific to particular theories but 
are somewhat universal in the sense that the same structure of divergences is common to 



QED, QCD and scalar QED in lowest order in the HTL resummation [2"I|- ||29|| . Further 



studies of the spectral function questioned the validity of the quasiparticle approximation and 
the exponential relaxation associated with a damping rate PU| , PT| . Although these studies 
provided an understanding of the failure of the quasiparticle picture (exponential relaxation) 
for hard fermions, the issue of the relaxation time scales was only recently clarified by the 
implementation of a Bloch-Nordsieck resummation of the infrared divergent diagrams |]32| , |33| 
which yields anomalous logarithmic relaxation. This resummation scheme was previously 
used at finite temperature to verify the cancellation of infrared divergences of soft photons 
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[34]. 

Infrared divergences in the propagators of charged fields are not particular of finite tem- 
perature field theory. It is well known that electrons in QED do not have a pole associated 
with their mass shell but rather a cut structure. This is a consequence of the emission of soft 
photons which because of their masslessness make the putative electron mass shell pole to 
become the beginning of a cut. This results in that even after ultraviolet renormalization, the 
wave function renormalization is logarithmically infrared divergent on shell. The implemen- 
tation of the Bloch-Nordsieck resummation of these infrared divergences at zero temperature 
leads to the correct electron propagator [35,36]. This resummation at zero temperature is 
equivalent to a renormalization group resummation of the leading infrared divergences in the 
Euclidean fermion propagator and leads to an anomalous scaling dimension (albeit gauge 
dependent) for the threshold behavior of fermions in QED [J35| . 

In this article we introduce a dynamical renormalization group resummation programme 
that allows to obtain the real time dependence of retarded propagators, and leads unequivo- 
cally to the real time dynamics of relaxation and thermalization without any assumptions on 
quasiparticle structure of collective excitations. This resummation scheme is the dynamical 
(real time) equivalent of the renormalization group resummation of divergences in Euclidean 
Green's functions which is so successful in both critical phenomena and asymptotic freedom 
and trascends approximations of the Bloch-Nordsieck type. The main concept in this pro- 
gramme is the resummation of secular terms in the perturbative solution of the equation 
of evolution of expectation values that determine the real time retarded propagators. This 
dynamical renormalization group was originally developed to improve the solutions of ordi- 
nary (and partial) differential equations |37| , and has been recently implemented in quantum 
field theory out of equilibrium |38| , |39[| where it reveals relaxation with anomalous (and non- 
perturbative) exponents (for other applications in quantum mechanics of few degrees of 
freedom see [f^Of] ). 

In the real time description of the dynamical evolution, the time variable acts as an 
infrared cutoff. The infrared singularities associated with the absorption and emission of 
massless quanta are manifest as logarithmic secular terms in the perturbative solution of the 
initial value problem. The dynamical renormalization group implements a non-perturbative 
resummation of these secular terms and leads to anomalous relaxation. In particular for 
scalar QED we find that the charged scalar field expectation value with hard momentum 
relaxes in absolute value as e - QTtln [*Ao] a ^ asymptotically long times. The asymptotic re- 
laxation is determined by the behaviour of the density of states p(k; u) as a function of u 
near threshold (u = \/k 2 + m 2 ). The larger is p(k;u) there, the faster is the decay of the 
expectation value of the field. 

The advantage of this method is that it leads to an understanding of relaxation directly 
in real time displaying clearly the contributions from different regions of the spectral density 
to the long time behavior. Furthermore, it offers a simple criterion to distinguish exponential 
relaxation and more complicated relaxational phenomena that cannot be interpreted within 
the quasiparticle picture. 

This method implements renormalization group resummations without the need for in- 
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voking a quasiparticle picture or any other approximation. 

The article is organized as follows: in section II we provide a direct link between linear 
response and the study of relaxation phenomena as an initial value problem out of equilib- 
rium. In section III we introduce and test the method of dynamical renormalization group 
within the simple setting of a field theory of two interacting scalars, one heavy and the other 
massless. This simpler theory presents the same type of infrared threshold singularities as 
scalar QED and QED. This model presents the infrared threshold divergences of a critical 
theory at the upper critical dimensionality. In this section we compare the Bloch-Nordsieck 
approximation and the renormalization group resummation of infrared divergences in the 
Euclidean propagator to the real time resummation implemented by the dynamical renor- 
malization approach in different situations at zero temperature. This study shows in detail 
the equivalence of all the different approaches at T = 0. We then implement the dynami- 
cal renormalization group at finite temperature and find anomalous logarithmic relaxation 
as in finite temperature QED [3"!2 ]. Section IV is devoted to a discussion of the dynamical 



renormalization group to elucidate this resummation program and to make contact with 
the usual renormalization in Euclidena space-time. In section V we study in detail scalar 
quantum electrodynamics. This theory has been previously studied within the imaginary 
time, equilibrium formulation and shown to have the same type of behavior as QED 
and QCD in leading order in the HTL resummation. We study both the exchange of bare 
photons and include the HTL resummation programme consistently to leading order into 
the dynamical renormalization group. This combined resummation of HTL and infrared 
secular terms in real time leads at once to anomalous logarithmic relaxation as in QED in 
the Bloch-Nordsieck approximation [32] and in the simpler scalar case studied in section III. 

We summarize our studies in the conclusion wherein we advocate to use this new approach 
based on the dynamical renormalization group to study fermionic excitations in a plasma 
and raise further questions and comments. The method of dynamical renormalization group 
leads directly to an understanding of damping and relaxation in real time without invoking 
a quasiparticle picture or any other approximation. 

Two appendices provide technical details and a third appendix provides a very simple 
and pedagogical example of the dynamical renormalization group. 



II. PRELIMINARIES: FROM LINEAR RESPONSE TO INITIAL VALUE 

PROBLEM 

We are interested in studying the real time evolution of expectation values of field op- 
erators. Consider a scalar field theory with an interacting Lagrangian density the 
expectation value of the scalar field $ can be obtained from linear response to an external c- 
number source term J. The appropriate formulation of real time, non-equilibrium dynamics 
is that of Schwinger-Keldysh [[I!]]- in which a path integral along a contour in imaginary 



time is required to generate all of the non-equilibrium Green's functions. 
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The non-equilibrium Lagrangian density along this contour is therefore given by |pH 
C NEQ [®+, J] = £[$+] + J$+ - £[$-] - J$- 

The non-equilibrium expectation value of the scalar field in a linear response analysis is 
given by 

/oo 
(Px*dt'G R {x -x,t- t')J(x, t') 
-oo 

with the retarded Green's function 

G R (x-x l ,t-t') = [G > (x-x > ,t-t') -G < (x-x l , t -t')]G(t-t') 
= i([<S>(x,t),<S>(x f ,t'])Q(t-t') 

where the expectation value is in the full interacting theory but with vanishing source. 
Consider an external source term that is adiabatically switched on in time from t — > — oo 
and of the form 

J(x',t') = J(x')e et 'e{-t r ) ; e^0+ (2.1) 

The retarded nature of G R (x — x',t — t') results in that 

(p(x,t = 0) = o (f) (2.2) 
<j>(x,t <0) = (2.3) 

where 4>o(x) is determined by J(x) (or viceversa, the initial value ^o(^) can be used to find 
J{x)) and the vanishing of the derivative for t < is a consequence of the retarded nature of 
Gr. The linear response problem with the initial conditions at t = given by Q2.2|) -( f2~3D can 



now be turned into an initial value problem for the equation of motion of the expectation 
value by using the (integro-) differential operator 0^,t) inverse of Gr(x — x',t — t') 

0^ t )(f)(x,t) = J(x,t) ; (f>(x,t = 0) = 4> (x) ; 4>(x,t < 0) = 



for the source term given by eq. ( |2.1| ). Within the non-equilibrium formulation the equation 
of motion of the expectation value is obtained via the tadpole method and automatically 
leads to a retarded initial value problem by coupling an external source that satisfies eq.(2~l). 



III. A SIMPLE EXAMPLE: A SCALAR THEORY 

We begin by considering a simple scalar theory of a massive and a massless scalar field 
with Lagrangian density 
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where the source coupled to the sigma field has been introduced to provide an initial value 
problem as explained in the previous section. Introducing the following renormalizations 



a 

2 rv 

m Z a 



6 J (T-r 



7T 



m 2 r Z a + 5m 2 ; g Z a Z x J 2 = g r Z g 



We now suppress the label r with all quantities being renormalized, and write the Lagrangian 
density in terms of renormalized quantities and counterterms 



C 



1 2 2 2 

—m a — qa ir 

2 y 



ct 



(Z a - 1) (<V) 2 - mV - -6m 2 a 2 + -iZ v - 1)(<V 



2 2 

-m a 



9(Zg 



l)a 2 TT 



--dmln 2 



where we introduced a mass counterterm for the 7r field to keep it massless. The counterterms 
are adjusted in perturbation theory as usual. 

The purpose of studying this simpler model is twofold: i) it provides a simpler setting to 
implement and test the method of the dynamical renormalization group and compare with 
previous studies of relaxation na t ure °f the infrared divergences in this simpler 

theory is similar to that of gauge theories in lowest order, i.e. the exchange of a massless 
field in the self-energy of a massive field. These infrared divergences are very similar to 
those found in scalar QED f29|, QED |2| and lowest order QCD 0. After studying the 
resummation of the infrared divergences via the dynamical renormalization group, we apply 
the method to gauge theories. 

Writing a(x,t) = if)(x,t) + (f>(x,t) with (ip(x,t)) = 0; 4>(x,t) = (a(x,t)) using the tad- 
pole condition and taking spatial Fourier transforms we find the equation of motion in the 
amplitude expansion 



4> k {t) + {Z a - 1) Mt) + u 2 k <j> k (t) + [u 2 k + 5m 2 } <p k (t) + / S fe (t - If) <p k (t') dt' = J k (t) 



with lo\ = k 2 +m 2 . We have absorbed the contribution of a momentum and time independent 
tadpole (ultraviolet and infrared divergent) in a renormalization of the mass. As discussed in 
the previous section, the source is chosen so that 4> k (t = 0) = <p k (0) ; cf> k (t < 0) = 0. £(£ — £') 
is the retarded self-energy. Writing Y> k {t — t') = -^rlk{t — t') (with the boundary condition 
7a,.(— oo) = corresponding to adiabatic switching-on of the interaction) and imposing that 
J k (t > 0) = 0, after an integration by parts the equation of motion for t > becomes 

4> k {t) + ( Z q - 1) \Mt) + u 2 k <f) k {t)] + \u>l + 5m 2 + 7fc (0)j <f) k {t) ~ f lk{t - t')Mt') dt' = 

L J L J _y 

(3.1) 



To one loop order we find 
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Sfc(t-t') 



471-3 



d 3 q 

q^ k+l 



{(1 + N q + sin [(g + a; fc+3 )(t - *')] 

(N q - n fe+9 ) sin [(g - u k+q )(t - t')}} 



(3.2) 



where N q is the Bose-Einstein distribution function for the massless field tt and n k+q is the 
corresponding distribution for the massive field a. It proves convenient to write the self 
energy in the form of a dispersion relation 



/oo 
du p(k; lo) sin [u(t — t')) 
-oo 

d 3 q 



(3.3) 



p(k;u) 
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47T 3 

consequently 



q 



{(1 + N q + n fc+(? )5 (u-q- u k+q ) - (N q - n k+q )5 (u — q + uJ k+q )} (3.4) 



Ikif) = I —p(k;u)cos(ujt) 



(3.5) 



A simple calculation yields 



p{k;u) 



2tt 2 



T , 
+ — In 



^ 2 

cu 2 — fc 2 
1 



+ — In- 



e t 



)(1 




_ <L_ 

e r 



(1 - e~V)(l - e 
[0(fc - w) + Q(k 



Q(u)- ui k ) 



(3.6) 



or 
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2(a; =f fc) 



There are some noteworthy features of the spectral density (|3.6|) above: i) whereas the zero 
temperature contribution vanishes (linearly) near threshold at uo = u k , the finite temperature 
contribution does not vanish at threshold, ii) the finite temperature contribution below the 
light cone —k<uo<k has its origin in Landau damping type processes in which the a 
particle scatters off a it particle in the medium. 

As it will be seen and understood below the spectral density ( |3.6|) leads to threshold 
infrared divergences. These will be studied in detail for the cases T = 0; T ^ separately 
in the next subsections. 

The equation of motion ( p.l| ) can be solved by Laplace transform. In terms of the Laplace 
transforms of <f) k (t) and 7fc(i) given by 4> k (s) ; j k (s) respectively, with s the Laplace transform 
variable, we find 



m 



s 2 + u 2 k + U{s) 



(3.7) 



with 
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C = l + (Z a -l)+[5m 2 + lk (0)] /ul 
U(s) = (s 2 + u\){Z a - 1) + 5m 2 + 7fc (0) - s %(s) 



(3.8) 



The term 



G(k,s) 



is recognized as the propagator in terms of the Laplace variable s. The retarded Green's 
function is obtained by the analytic continuation 

G ret (k, to) = G(k, s = iu + e) | e=0 + 

The Laplace transform of the self-energy is recognized to be 

Sfc(s) = 7fc(0) - s 7 fc (s) = J dup{ 

and its analytic continuation is then given by 

t(s = ioo + + ) = E R (u) + i£/(w) 



LU 



S 2 + LU 2 



Z r (uj) = J duj'p{k,uj') V 



L)' 2 — id 2 



TV 



(3.9) 
(3.10) 



therefore LT (s) is recognized as the twice subtracted self-energy which is rendered finite by a 
proper choice of counterterms. Furthermore, choosing to renormalize at s 2 = —uu 2 with the 
counterterms given by 



Z a -\ 



d^njk , to) 
dco 2 



\U)=Ul k 



5m 2 = - [ p(k,u') V " 2 dJ 
J uo 2 - uoi 



we find 



1 + 



/ 



LO 1 



V- 



J 2 - u 2 ) 2 



dJ 



(3.11) 



is finite even for renormalizable theories in which p(k, to) ~ lo 2 at large to leading to quadratic 
and logarithmic divergences with logarithmically divergent wave function renormalizations. 
This finite wave function renormalization will be seen to emerge naturally from the dynamical 
renormalization group. 

The real time evolution is obtained by performing the inverse Laplace transform along a 
path in the complex-s plane parallel to the imaginary axis to the right of all the singularities 
of 4>k(s). We note that the putative pole at s = has vanishing residue. 
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A. Bloch-Nordsieck and Renormalization Group at T = 



At T = we find after renormalizing the mass 



= 7 fc (0) - S7fc(s) 



g 2 Pl + m 2 1 P% + m 2 



P 2 

E 



4tt 2 
s 2 + k 2 



PI 



(3.12) 



where P§ is identified with the Euclidean four momentum squared. Since in this theory the 
wave function renormalization is finite, we choose Z G = 1 in what follows. Therefore up to 
this one loop order we obtain the Euclidean irreducible two point function (the inverse of 
the Green's function) to be given by 



T(k,s) = [G(k,s)y 




1 + 



In 



p2 



+ 1 



m 



R 



(3.13) 



We clearly see that P\ — > —m 2 R is no longer a pole but the end point of a logarithmic branch 
cut, corresponding to the threshold for the intermediate state of a a particle and a soft 
massless 7r particle. 

This logarithmic infrared divergence near threshold is the same phenomenon as in gauge 
theories wherein the self energy of charged fields has an infrared logarithmic divergence at 
threshold associated with the emission of soft quanta. The Bloch-Nordsieck resummation 
exponentiates the logarithmic divergences near threshold leading to 



G BN (k,s) 



1 



P 2 



A-l 



+ 1 



m 



R 



m 



(3.14) 



with the dimensionless coupling 



A 



9 Z 



^K 2 W? R 



(3.15) 



Using this resummed expression for the propagator the inverse Laplace transform can be 
performed by wrapping the contour around the branch cuts along the imaginary axis from 
s = ±iuk to ±ioo. Computing the discontinuity of GBN^k, s) across these cuts the real time 
evolution of the expectation value is given by 



0fe(O) 



LulCm 2 ^ -sin[n\] [ 

7T Ju: 



00 du cos(uot) 



21 1-A 



with io\ 



k 2 + m 2 . 
Using now the result [[46j 

dx cos(xy) 



\1-X 



k 2 y 
4 sin7rA 



{ji_ A (y) [N_i_ x (y) 



H 



9 



J^M [N hx (y)-m__ x (y)]} 



with J v {z) a Bessel function and TC u (z) a Struve function, we find the asymptotic long time 
behavior 



m 2 r(l-A) 



m 

4 



l-A 



2 / vrA N 

cos I Uk t H — 



1 + 



/ 1 



\UJ k t 



(3.16) 



Here we used the asymptotic formula 
H v {z) - N v {z) z ~-- 









(r 


1 - 


"Gil 









We see that the Bloch-Nordsieck resummation of the infrared divergences leads to relaxation 
with an anomalous exponent. This is similar to the case of QED. We now argue that the 
infrared divergence and the emergence of anomalous dimensions can be understood by es- 
tablishing a parallel with static critical phenomena at the upper critical dimensionality in 
Euclidean space-time. This connection will pave the way to using the renormalization group 
to sum up infrared divergences non-perturbatively much in the same manner as in the theory 
of critical phenomena. In order to establish this connection more clearly we now introduce 
the dimensionless variable 

^2 Pi 



P 



m R 



the main reason for introducing this variable is that when P| — > —m 2 R ; P ^0, therefore 
the threshold behavior is mapped onto the zero Euclidean four momentum region in terms of 
the new variable. In critical phenomena logarithmic divergences appear when the Euclidean 
four momentum goes to zero at criticality at the upper critical dimension. We now introduce 
a wave function renormalization constant 

Z^(K) = 1 + WnK 2 

and a renormalized irreducible two-point function 



m 2 R P 2 



1 + A 




lnP 2 + lnK 2 



T R (P,K)=Z <t> (K) T(P) 

where A is the dimensionless coupling ( |3.15| ) and K is an arbitrary renormalization point. 
The bare irreducible function T(P) is independent of the renormalization point K , i.e. 
KdT/dK = which leads to the renormalization group equation 



= dZ^K) 
<91nZ 



T R (K,P) = 



2A 



(3.17) 



10 



Near threshold when P — > 



Tr{K,P) 



m, 



P z 




with $ a dimensionless function of its argument. The renormalization group equation ( |3.17 ) 
then leads to 



with solution 



~n 9 

p w +r] 








p 
T 



finally leading to the renormalization group improved two-point function 



G RG (k,s)= ®{l)m\K 



-i 



p2 



1-A 



(3.18) 



which coincides with the one obtained by the Bloch-Nordsieck resummation, eq. ( 3.14Q up 
to an overall multiplicative factor. We can now retrace the same steps that led to the real- 
time evolution of the expectation value, by performing the analytic continuation for the 
retarded correlation function and the inverse Laplace transform, leading to the relaxation 
with anomalous dimension given by eq.( |3.16|) . 

The equivalence between the renormalization group improved and the Bloch-Nordsieck 
re-summed propagator is fairly well known |35|j36|| . The main purpose of our analysis is to 
make the point that the anomalous dimension in the amplitude of the expectation value can 
be understood as arising from the scaling behavior of the Green's function near threshold in 
terms of the variable P. This scaling behavior, a result of the infrared divergences associated 
with the emission of soft quanta is akin to those in static critical phenomena. 

Furthermore, this structure of anomalous dimension of the Euclidean propagator as a 
result of threshold infrared divergences is similar to that found in QED at T = either via 



a Bloch-Nordsieck or renormalization group resummation [35,36]. 

In preparation to the forthcoming discussion presented below on the real time interpreta- 
tion of the renormalization group resummation of infrared divergences, it proves illuminating 
to obtain the perturbative form of the real time solution. This will allow us to identify the real 
time manifestation of infrared divergences. The naive perturbative expansion of the renor- 
malization group improved propagator ( |3.14| ) to first order in the coupling A near threshold 
leads obviously to the one- loop result ( |3.13| ). After performing the Fourier transform and 
obtaining the real time evolution given by (|3.16|) we can now expand naively in the coupling 
constant, and find 



4> k (t) « A k cos (u k t) [1-A ln(uj k t)] + 



(3.19) 
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with A k the amplitude read off from ( |3.16| ). This expression reveals that if we attempt a 
perturbative solution of the real-time equation of motion ( |3.1|) we would find logarithmic 
secular terms, i.e. terms that grow in time and invalidate the perturbative solution at long 
times. In this case a perturbative expansion of the real-time equation of motion will break 
down at time scales tb reak ~ e 1//A /uj k . The renormalization group in the energy representation 
provides a resummation of the infrared divergences in the propagator, which leads to a real- 
time evolution that is asymptotically decreasing function of time. 

We now study the perturbative solution of ( |3.1| ) that reveals indeed these secular terms, 
and implement a real-time version of the renormalization group that implements precisely 
this resummation. 



B. Dynamical Renormalization Group at T = 

Having established the resummation of threshold infrared divergences both within the 
Bloch-Nordsieck approximation and the renormalization group, we now introduce a novel 
method that allows a similar resummation but directly in real time. Consider seeking a 
solution of the equation of motion ( |3.1| ) in perturbation theory in the coupling. Writing 
the self energy as an expansion in terms of the dimensionless coupling A given by ( 3.15| ), 



^fc = J2 n =i a perturbative solution obtained as a power series expansion is given by 

4>k{t) = 4> k °\t) + (0 + " " ' with the hierarchy of equations 

4> k 0) (t)+^ k °\t) = o 

4 13 W + (* ) = - [8m 2 + 7? } (0)] 0f (t) + f l { k 1] (t - 4 0) (0 dt< 

1 J Jo 



where we note that the contribution from the wave-function renormalization vanishes by 
virtue of the zeroth-order equation of motion. The solution to the zeroth-order equation is 

0<W(t) = A k e iuJkt + A* k e~ iuJkt 

and the initial conditions <f) k (t = 0) = 0&(O); ;<f) k (t = 0) = implies A k = A* k = <ft k (0)/2, 
but we will leave both constants to recognize more easily the different contributions and we 
will use this condition at the end of the calculations. The solution to the above hierarchy of 
equations can be found in terms of the retarded Green's function of the unperturbed problem 

Qii(ti - t 2 ) = — sin[u; fc (ti - £ 2 )]©(*i - h) 
u k 

The higher order corrections are easily (but tediously) computed using the spectral repre- 
sentation of the self-energy. The first order term is given by 
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^' a \t) = — ( —p(uj) f ' dt' C dtx sin[u k (t-t')]cos[uj(t' -tx) (o) (*i; 

LU k J LO JO Jo 



— (dm 2 + / 
u k \ J 



—p{uu)\ t smluJt - t')U (0 \t') dt' 
to Jo 



where the spectral density p(to) is given to one loop order by eq. ( |3.4j ) and its leading infrared 
contribution for T 7^ given by eq. (|3.6| ). The integral over the time variables can be done 
straightforwardly, the result is given by: 



A J lu 

A k e~ iu}kt 



p(u) <^A k e 
1 



t — I +UJ 



UJ k - LU \ %{U- LU k ) 

e 2ico k t _ ^ e i(ui k +w)t _ ^ N 



-LO 



LU k - LU 



A* k e iUkt 



A* k e- iuJkt 



1 



2iu k 

-2iu k t 



i(uu k + to) 



+ LO 



-lu 



+ 



LU k + LO \ 
1 

UJ k + to 



-H101 



\ e ~i{u k -u)t _ ^~ 
- + 



i(u k - uj) 



LU — > —LU 



' e i(u+w k )t _ ]\ 

t ~ — +LU 

l(LU + LU k ) I 



-LO 



(3.20) 



i ( . 2 r duu 
— (5m 2 + — 

2uu k \ J LU 



P (lu) [A k e %u>kt t - A* k e~ iu>kt t 



A k e~ luJkt 



„2iw k t _ 1 -2icj k t 

_ £ _ A * ^ k t e 



1%LUi 



IlLUl 



(3.21) 



Secular terms will arise from the contributions of the form e 4 ^""^)* — 1 if the coefficients of 
these terms produce singularities in the integration region. 

We are now in condition to analyze different cases. Although we are primarily interested 
in applying the method of dynamical renormalization group to the situation of infrared 
divergences, in order to gain insight and test this method we begin by studying situations in 
which results are known. To this purpose we address the familiar case of a generic interacting 
scalar theory in which the pole frequency to k is away from thresholds {outh}-, either above, in 
which case there is a resonance, or below in which case the particle is stable. 



1. LO k ^ LUth 

For u k < Lu t h there are no singularities in the integration region, therefore the only secular 
terms are those linear in time in ^'^(t) ; (j) k (t) . If on the other hand uu k > LU t h and far 
away from threshold, there are singularities (simple and double poles) in the integration 
region. 
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We can extract the secular terms in the above expression in both cases uj k > Wth an d 
Uk < uJth by using the results of appendix B or alternatively taking the long time limit using 
the distributions (see appendix B) 



lim 



t — > oo — [at — sin at] = V I — ] [at — sin at] 



of 

lim 1 

t — > oo — 



1 



a< 



cos at 



a- 



wtS(a) + V 



a' 



(1 — cos at) 



where the 5(a) accounts for resonant denominators. This term is recognized from the familiar 
Fermi's Golden rule. 



Gathering the secular terms from both contributions ( |3.20|J3.2l| ) and taking A k = A* k we 



find 



^\t)=A k e^< 



,8m 2 + E R (uJk)_ _ S/MQ 



2uj k 



+ 



duj 



p(u) V- 



[OJt 



10 



2\2 



+ c.c 



(3.22) 



where we have used the expressions for the real and imaginary parts of the analytically 
continued self-energy given by equations (|3.9|) - (|3*n0 ). The last terms are non-secular at long 
times and remain perturbatively small. 

In this manner the resummation of the secular terms is obvious and correspond to a shift 
in the pole position Uk — * Uk+5uJk (finite by a proper choice of 8m 2 ) and a width or decay rate 
Yk given by the imaginary and real part respectively of the first order correction above. Since 
this is the simplest and most familiar setting to introduce the dynamical renormalization 
group, we now present the resummation of the secular terms via this method. This is achieved 
by introducing a (complex) renormalization of the amplitude and writing P^JSH 



A k = A k (r) Z k (r) 
Z k (r) = l + \z^{T)+i\zf\r) + 



(3.23) 
(3.24) 



where r is an arbitrary time scale that acts as a renormalization point and the Zr\(t) are 
real functions. Choosing 



A*£>(r) 



£j(^fc) 



t; Xz?\t) 



8m 2 + T> R (uj k ) 
2luu 



we obtain 



.8m 2 + H R (u k ) E/(w fc ) 
1 + 1 ~ \t — t) n \t 



2ui 



2ui 



T 



+ c.c + regular terms 

(3.25) 
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where the regular terms refer to the non-secular last term The meaning of the above expres- 
sion is clear: a change in the time scale corresponds to a change in the (complex) amplitude 
of the expectation value. Whereas the original perturbative expansion was only valid for 
times such that the contribution from the secular terms remain very small compared to the 

remains valid for intervals t 



unperturbed value, the renormalized expression ( |3.25| ) remains valid for intervals t — r such 
that the secular terms remain small. By choosing r arbitrarily close to t we have improved 
the perturbative expansion. However 4>k{t) does not depend on r: a change of the renormal- 
ization point r is compensated by a change in the complex amplitude A k - This leads to the 
dynamical renormalization group equation to lowest order 



dAkjr) 
dr 



.5m 2 + E R (oo k ) X/Ofc) 



2ujk 



A, 







with obvious solution 



A k (r) 
5oo k 



A k (0) 
5m 2 + T, R (u k ) 
2oo k 

2u k 



Choosing the arbitrary scale r to coincide with the time t we obtain the resummed 
expression for the expectation value 



<f> k (t) = CA k (0) e 



C = 1 + 



duo 



■ -< l -A k {()) f—e 
J oo 

ool 



p{u) V- 



"1 



ul-u 2 ) 2 



+ c.c 



(3.26) 



CO 



p{uo)V 



with ujp(k) the pole position shifted by one-loop corrections and T k is identified with the decay 
or damping rate. The constant C is the same as in eq. ( 3.11 ). The on-shell renormalization 
leading to eq. ( p.ll[ ) is here a consequence of the perturbative expansion in terms of the 
solutions of the equations of motion. 

After some straightforward algebra, the constant C is found to be the same as the residue 

at the pole (or resonance) at uo p . The last, non-secular terms 



of the Laplace transform 
in 



3726], allows us to make contact with previous results j|4j . The long time dynamics of this 

If the spectral density vanishes near 

oo 



integral is dominated by the threshold contribution |44 
threshold as p(uo 
law relaxation t~ 



ooth) a then the asymptotic time evolution is described by a power 
1 (long time tails). Thus we see that the dynamical renormalization group 



resummation has obtained all of the features of the solution via the Laplace transform (|3.7| ) 
which were previously obtained jPJ. 

The resummation via the dynamical renormalization group has led to a (asymptotic) 
convergent perturbative expansion for the time evolution of the expectation value. 

This simple case provides for a clear understanding of the resummations implied by the 
dynamical renormalization group and paves the way for understanding the more complicated 
cases of threshold singularities and finite temperature below. 
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2. Threshold singularities 



Having established the reliability of the dynamical renormalization program in more 
familiar settings, we are now in position to apply this method to study the case of threshold 
infrared divergences arising from the emission of soft massless quanta. Thus we now return to 
the theory of a massive and a massless scalar fields of the beginning of this section. We begin 
by analyzing the situation at T = to make contact with the Bloch-Nordsieck and Euclidean 
renormalization group resummations, but now implementing the dynamical renormalization 
group resummation. 

Since near threshold the spectral density 

(El)) 

p(k, uo) — k — 4 A oo k {oj 



for T = becomes (A is defined by eq. 



wjb) + O{uo 



(1,6) 
k 



besides the linear secular terms in 
tentially infrared divergent term arising from 
appendix B, found to be given by 



(see equations ( p.2U| )-( |5^TD there is a po- 
which is real and using the results of 



duo 



CO 



p(k; uo 



1 - cos[(u; - U k )t] M t>>] 
(a; - oo k ) 2 



-A(ln/it + 7 ) +0{\) 



where 7 is Euler-Mascheroni constant (see details in appendix B). We note that time acts as 
an infrared cutoff in the sense that for 00 ~ co k at finite time the integral is convergent. The 
infrared divergences are now manifest in a logarithmic time dependence. The linear secular 
terms combine just as in the previous case to provide an imaginary secular term given by 



i(5m + T lR (to k ))/2io k just as in eq. (|3.22| ) which in this case is simply frequency shift as can 



be seen from the expression for the self-energy given by eq. (|3.12|) . This shift is made finite 
with a proper choice of Sm 2 . Thus in this case we find 



.(*) 



A k e luJkt 



1+i 



Sm 2 + T, R (u k ) 
2io k 



t — A In jit 



+ c.c + regular terms 



with p = /ie 7 . Similarly to the previous case, we introduce the complex amplitude ( p. 24 ) 
and choose 



a4 1} (t) = Aln/Zr ; AzfV) 
leading to the following expression for 4>k{t) 



Sm 2 + T, R (u k ) 
2oo k 



-S00 k T 



Mt) = Mr) e^ kt 



1 + i5co k (t — r) — A In — 



+ c.c + regular terms 



4>k{t) is independent of the arbitrary time scale r leading to a renormalization group equation 
obeyed by the complex amplitude, which is now given this order by 



dA k {- 
Or 



iduoi 



A" 



T 



Mr) = 



16 



with solution 



A k {r) = A k (r ) e 



To 

T 



i A 



Again, choosing the scale r to coincide with the time t, we finally obtain the asymptotic 
dynamics of the expectation value in this case to be given by 



TO 
t 



+ c.c + small 



(3.27) 



where u>k,R — + 5u>k- The terms denoted by small remain perturbative at all times and 
decay faster than the term with the anomalous dimension in weak coupling. This expression 
coincides with the long time behavior found in the previous sections via the Bloch-Nordsieck 
and the renormalization group resummation of the logarithmic infrared divergences of the 
propagator given by eq. ( |3.16| ). We thus conclude that the dynamical renormalization group 
implements a resummation in real time which is complementary to the renormalization group 
or Bloch-Nordsieck resummations in the frequency representation of the propagator. 



C. Dynamical Renormalization Group resummation 



At finite temperature the infrared divergences are enhanced by the Bose-Einstein distri- 
bution function of massless particles N q « T/q ; q/T « 1. This can be seen at the level 
of the spectral density p{k;u) given by eq. (|3.6p . Whereas the zero temperature contribution 
vanishes linearly near threshold, the finite temperature contribution remains constant there 
provided m ^ 0. In particular we find that near threshold the Laplace transform of the 
retarded propagator behaves as 



T{k,s;T) 

g(T,k) = 



2 i 2 
E + ™R 



1 + 



g 1 



4tt 2 P! 



In 



Pi 



+ 1 



m 



R 



-g(T,k) In 



Pe + ™ 2 r 
/i 2 



4tt 2 



In 



+ k 

LUfc — k 



where \x is an arbitrary infrared cutoff scale and assumed that the external momenta is 
not too hard. Thus whereas the zero temperature inverse propagator actually vanishes 
at threshold (with an infrared divergent slope) the finite temperature propagator diverges 
there, reflecting the stronger infrared divergence at finite temperature. In this situation a 
Bloch-Nordsieck resummation of the Euclidean propagator is not clear as was emphasized in 



reference [32| and a multiplicative (wave function) renormalization cannot cure the infrared 
divergence since the finite temperature part is not proportional to P| + m R . 

It is precisely in this situation that the power of the dynamical renormalization group is 
revealed. We will study two different cases in detail: i) hard external momentum k » m 
(or m — 0) and ii) soft external momentum k < in. 
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1. Hard external momentum (m ~ 0) 



In this case the Landau damping cut coalesces with the cut for u > k and both terms 
of the spectral density ( |3.6| ) contribute. Because of the simmetry uj — > —uj of the time 
dependent terms in eq.( |3.20|) the frequency integral in the interval —k < uj < k plus the 
integral in the interval k < uj < oo can be folded into an integral in the range < uj < oo in 
terms of the effective finite temperature spectral density 



p(k;u) 



l_T_ 

TT 2 k 



In 



2T 



en 



The asymptotic time dependence is dominated by the region uj « k which gives the infrared 
divergences of the propagator. In this region the effective spectral density 



p(k; uj) "= k In e(o,) + 0(uj - k) 



T 



_k_ 

e t 



(3.28) 



We start by analyzing the different contributions to the infrared behavior of the coefficient 
of A^e tulkt . A simple analysis leads to the following conclusions: i) the contribution near 
threshold to the imaginary part of the coefficient cancels out between the production cut 
(uj > k) and the Landau-damping cut (0 < uj < k) leaving a linear secular term without 
infrared divergences that renormalizes the mass, ii) The contribution near threshold to the 
real part is the same for both cuts < uj < k and uj > k and add up. This contribution in 
the asymptotic long time limit is obtained from the formulae in appendix B and given by 



T 



n 2 2k 2 



t 



f 00 dz . 


z 


/ ~2 ln 




/o Z A 


.2tT. 



1 



COS z 



Aixk< 



T t In/It ; JZ = 2T e 7_1 



(3.29) 



where we have quoted the leading contribution in the asymptotic time regime. Subleading 
terms can be consistently obtained using the formulae of appendix B. The coefficient of 
A* k e lu)kt has a potential infrared divergence, however the contribution from both cuts cancel 
each other leaving an infrared (and ultraviolet) finite result without secular terms. The zero 
tempreature contribution and that of (^'^ lead to a finite frequency shift. Implementing the 
dynamical renormalization group resummation we find 



(f) k {t) = A k {t ) e^* ( '-' o) e~ akTt ln[t/ta] + c.c. ; 



4nk 2 



where we have solved the renormalization group equation with an initial condition at a time 
to = This solution reveals clearly the renormalization group invariance, a change in the 
arbitrary time to is compensated for by a change in the amplitude and an overall phase. 
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2. Soft external momentum (m ^ 0) 



For m^fl the Landau damping cut and the production cut are separated and the infrared 
divergences arise only from the production cut uj k < uj. In the high temperature limit we 
find near threshold 



p{k;uj] 



7T" 



T 
~k 



\nM k [1 + O{uo - u k )} Q(uj - uj k ) 



M k = 1 



2k n \ 
— (k + UJ k ) . 



ni 



(3.30) 



We note that in this case the spectral density p(k; to) approaches a constant value at 
threshold. 

The terms proportional to t/(u + uj k ) in ^k' a \t) do not have infrared divergences but 
they remain as secular terms that combine with those of (t) to give a renormalization of 
the frequency just as in the previous cases. The infrared divergences arise from terms with 
denominators l/{yj — uj k ). 

In this case these infrared divergences are manifest as logarithmic secular terms in the 
real and imaginary parts leading to damping and anomalous logarithmic phases. 

These contributions are the following: i) the imaginary part of the coefficient of Ae tuJkt is 
given asymptotically by (see appendix B) 



i t 
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+i t 



uik UJ UJ — Ufa 



1 



duo 



uj k UJ — UJk 



p{k)u) 



sin(u; — Uk)t 
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p(k\uj k ) 



nt»i ig 2 lnM k 



UJ 



Q(p — uj + uj k 




T t In [/it e 



7-ii 



(3.31) 



Notice that eq. (|3.31 ) is independent of the scale p as one can easily see since the derivative 
with respect to p of the r. h. s. identically vanishes. The scale p has been introduced just to 
have a dimensionless argument in the logarithms, ii) The real part is asymptotically given 
by (see appendix B for details) 
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,2 



COS w 



UJ 



UJ k )t] 



1 



1 



9*Tt 
8 7T kuj k 



9 2 T 
A7T 2 m 2 uj\ 



1 + 



{UJ - UJ k 

3 m 2 



fit»i 



2kuj k 



\nMi 



2u k (uj - uj k )_ 
\n\pte^] + 0{g 2 ) 



The remaining secular but infrared safe terms in 0( 1,a ); contribute to the imaginary 
part a term proportional to t which can be absorbed in a redefinition of the arbitrary scale 
p in eq. (|3.31|) . The logarithmic divergence as k/m —>■ oo reflects precisely the logarithmic 
time dependence found in the previous case of hard momentum m ~ and that results in 
the anomalous relaxation proportional to t In pt as given by eq. ( |3.29| ) . 

Implementing a resummation of the secular terms with the dynamical renormalization 
group following the steps outlined above leads to the asymptotic form 
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+ c.c. 



5/c 



£ 2 T 
47r 2 m 2 cl^ 



InM 



3m2 i a/ 
1 + 771 lnM * 

ZKUJk 



(3.32) 



It is interesting to try and understand the exponential damping in this case. For this we 
show in figure 1 the spectral density for the sigma field in this soft case near threshold 



S(k, u; T) 

E R (u, k; T) 
S/(a;, k] T) 



E T (uj,k;T) 



E R {cu, k; T)Y +[Ej(w, A:; T)\ 



-p(k,u k ) In 



+ O{uo-uo k ) 



p{k,u k ) + 0(u - LO k ) 



(3.33) 



where p(k,u k ) is given by eq.( p.30[ ) and for the figure it was taken to be p(k,u k ) = —0.005. 
S(k,uj;T) vanishes at threshold but with a singular slope, this results in that the spec- 
tral density features a sharp peak near threshold, which is found to be at uo w uo k + 
\p{k,u k ) In \p(k,u k ) \ H . 

The decay rate is given by (l/2)E/(o;/ c )/(2a;fc) the extra factor 1/2 is simply a result of 
the fact that u k is the threshold, i.e. the end-point of the integral and therefore the on-shell 
delta function only picks-up half of the contribution. The logarithmic dependence of the 
phase is a consequence of the logarithmic infrared threshold divergences and prevents an 
interpretation in terms of quasiparticle poles. 

This case must be contrasted to that of m = 0, wherein the imaginary part of the self- 
energy is infrared singular at threshold. This is revealed in the logarithmic singularity in 
the limit k/m —>■ oo which reflects precisely the logarithmic time dependence found in the 
previous case of hard momentum m ~ leading to the anomalous relaxation proportional 
to t In [pit] as given by eq. ( |3.29| ). 



IV. DISCUSSION AND INTERPRETATION: 

Before proceeding further to the case of a gauge theory, it is convenient to pause and 
analyze the results that we have obtained so far and elucidate the main aspects of the dynam- 
ical renormalization group. The scalar model chosen in the previous section is a non-trivial 
example of a superrenormalizable theory that displays the same type of infrared divergences 
as a critical theory at the upper critical dimension. The usual renormalization group leads 
to a Bloch-Nordsieck resummation and anomalous dimensions at T = 0. Performing the 
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Fourier transform we recognized that the real-time interpretation of the renormalization 
group resummation corresponds to a power law relaxation with the power determined by 
the anomalous dimension. A naive perturbative expansion in the dimensionless coupling 
results in secular terms, i.e. terms that grow in time and signal the breakdown of the pertur- 
bative expansion at long time scales. These same secular divergences are obtained directly 
in real time when the equation of motion for the expectation value is solved in a perturbative 
expansion. The dynamical renormalization group implements a resummation of these secular 
divergences which leads to an improved perturbative solution. The renormalization proce- 
dure can be understood with a very simple and pedagogical example, the weakly damped 
harmonic oscillator. Consider the equation of motion 

V + y = -ey , e << 1 

attempting to solve this equation in a perturbative expansion in e leads to the lowest order 
solution (see appendix C) 



where the term that grows in time, i.e. the linear secular term leads to the breakdown 
of the perturbative expansion at time scales tb rea k oc \je. The dynamical renormalization 
introduces a time scale r in the form A = A(t) Z{t) ; Z{t) = 1 + e z±(t) + ■ • •, choosing z\ 
to cancel the secular term at this time scale leads to the renormalization group equation 



and the improved solution y{t) = e~^ t (A(0)e lt + c.c) after setting r = t in the solution. This 
obviously is the correct solution to 0(e). The interpretation of the renormalization group 
resummation is very clear in this simple example: the perturbative expansion is carried 
out to a time scale r << 1/e within which perturbation theory is valid. The correction is 
recognized as a change in the amplitude, so at this time scale the correction is absorbed 
in a renormalization of the amplitude and the perturbative expansion is carried out to a 
longer time but in terms of the amplitude at the renormalization scale. The dynamical 
renormalization group equation is the differential form of this procedure of evolving in time, 
absorbing the corrections into the amplitude (and phases) and continuing the evolution in 
terms of the renormalized amplitudes and phases. This is the same spirit as the momentum- 
shell renormalization in critical phenomena. The details of the second order calculation and 
implementation of the renormalization group for this simple problem are offered in appendix 
III to illustrate the shift in the frequency. 

This interpretation in the simple exercise extends to the more complex situations with the 
same underlying mechanism, the secular divergences are absorbed in the complex amplitudes 
and the perturbative expansion is then carried in terms of the renormalized amplitudes. The 
differential form of this process is the dynamical renormalization group equation. Thus the 
similarities with the usual renormalization procedure are manifest. 



y(t) = A e il 1 



--tj+c.c 
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In the same manner that the usual renormalization group sums the leading logarithms 
when the renormalization group functions are computed to one loop order, the dynamical 
renormalization group sums the leading secular terms when the coefficients are computed 
to lowest order. This is manifestly revealed by the naive perturbative expansion of the real 
time propagator ( |3.16| ). 

The main objective of studying the different cases in the scalar theory of the previous 
section was to thoroughly test the method with a non-trivial example that can be studied 
with different methods. By studying these examples in detail we have learned that the dy- 
namical renormalization group provides: i) a real-time equivalent of the resummation via the 
renormalization group in Euclidean space time (equivalent to Bloch-Nordsieck resummation) 
in the case of infrared threshold divergences at T = leading to relaxation with anomalous 
exponents, ii) the usual mass shift and damping rates in the case of narrow resonances and 
iii) leads to a similar resummation scheme at T ^ when the threshold infrared divergences 
are more severe. This detailed analysis then provides confidence on this novel method to 
study the more interesting and relevant case of a gauge theory. 



V. A GAUGE THEORY: SQED 

We are now in position to apply the method of the dynamical renormalization group to 
implement the resummation of infrared divergences in gauge theories, which is our primary 
goal. 

We will study the case of scalar QED, since to lowest order in hard thermal loops this 
theory has the same properties as those of QED and QCD [(10,18,^9,321, in particular the 



infrared divergences associated with the propagation of the charged fields. 

Since we are primarily interested in studying the real-time manifestation of the finite 
temperature infrared divergences, we will focus on the relaxation of the charged scalar field 
at finite temperature. Furthermore we will only consider the contribution of transverse 
photons to the charged scalar self energy, since longitudinal photons are Debye screened at 
finite temperature (m^ « eT) and do not contribute to the infrared divergences. 

In this Abelian theory it is rather straightforward to implement a gauge invariant formu- 
lation by projecting the Hilbert space on states annihilated by Gauss' law. Gauge invariant 
operators can be constructed and the Hamiltonian and Lagrangian can be written in terms 



of these. The resulting Lagrangian is exactly the same as that in Coulomb gauge |^5[ and is 
given by 

£ = fy$t - m 2 &<3> + -O^At ■ d»A T - eA T ■ j T - e 2 A T ■ A T $ f $ + 
+ ^(VA ) 2 + e 2 A 2 &$ + eA p, 



j T = i($tv r $ - Vr^) ; p=-i ($$ f - $ t $) . 



where we have traded the instantaneous Coulomb interaction for a gauge invariant Lagrange 
multiplier field A which should not be confused with a time component of the gauge field. 
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At is the transverse component satisfying V • At(x, t) = 0. Since we are only interested in 
obtaining the infrared behavior arising from finite temperature effects we do not introduce 
the renormalization counterterms to facilitate the study. The finite temperature behavior 
is ultraviolet finite. The non-equilibrium generating functional requires the fields on the 
forward and backward branches f45fl . The equation of motion for the charged scalar field is 
obtained by writing 



$ ± (x,t) = <p(x,t) + A ± (x,t) ; (A ± (x,t)) = 

and similarly for the hermitian congujate fields. 

In obtaining the equation of motion to one-loop order, we neglect the contribution from 
the Coulomb interaction. The reason as explained above is that the long-range Coulomb 
interaction will be screened by finite temperature effects with a Debye screening length 
mo ~ eT and hence the screened Coulomb interaction will be free of infrared divergences. 
However in an Abelian plasma the magnetic (transverse) photons are not screened (no static 
screening, only dynamical screening through Landau damping) and the exchange of soft 
magnetic photons will lead to threshold infrared divergences. 

In terms of the spatial Fourier transform of <p(x, t) we find the equation of motion 



<p{-k, t) + (k 2 + m 2 + e 2 (A 2 T (x , t)))<p{-k, t) - Ate 2 jH dt' J -0- k T {q) k T (q) x { 

(Ap(q, t)Ap(-q, t'))(A+(-k - q, t)A^+(k + q, t'))- 
(Ap(q,t)A^(-q,t>))(A + (-k - q,t)A^-(k + q,t>))} <p(-k,t')} = J{-k,t) 
k T {q) = k- q(k-q) 

where we coupled an external source J(k,t) = J(k)e €t Q(—t) (e — > + ) to provide an 
initial value problem with an adiabatic switching on of the expectation value. The initial 
conditions for the expectation value are 

(p$,t = Q) = ¥>*(0) ; <p(k,t = 0) = 

Infrared phenomena is associated with the soft limit of the intermediate photon, and 
therefore requires the HTL resummation of the intermediate photon propagator. However, 
we begin our study of SQED by implementing the dynamical renormalization group to 
resum the perturbative expansion in the case in which the self-energy of the charged field 
only includes the exchange of a bare transverse photon. The study of this situation will 
shed light on the different physical processes that contribute and those that do not because 
of dynamical screening via Landau damping. In the next section we will include the HTL 
resummation of the exchanged photon and implement the dynamical renormalization group. 

A. Bare photon propagators 

The necessary free-field non-equilibrium Green's functions are given by 
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•Scalar Propagators 
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G k -(t,t') = G>(t,t')e(t' -t) + G<(t,t')Q(t - t') , 
G^(t,t') = -G< > \t,t'), 
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•Photon Propagators 
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Here Vijik) is the transverse projection operator: 



(5.1) 
(5.2) 



Finally we find the equation of motion for t > to be given by 

£(M -0 <p(-k,lf) dt' = 

-oo 

E (*, t - f ) = -8e 2 k 2 J ,f«\~™ 26 {(l + N q + n u .) sin(^ +(? + q) (t - t>) 



W 4 Q 
+ {N q - n Uq ) sin(^- q)(t - t')} 
M 2 (T) = m 2 + e 2 < A 2 > 
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and 9 is the angle between k and q. We see that the self-energy has a form very similar to 
that of the scalar case, eq. ( |3.2|) , the only difference being the k 2 in front (reflecting the 
exchange of transverse photons) and the 1 — cos 2 9 inside the integral. Integrating by parts 
and using the initial conditions the resulting equation of motion can be written in the same 
form as in eq. (|3.1|) with 7fc(t) as in eq. (|3.5|) in terms of the spectral density 



e 2 k 2 



p{k;u) 



2tt 2 



qdq 

^k+q 



1 — cos 2 9) dcos9 



{(1 + N g + n k+q ) 8(u-q- u k+q ) 
{N q - n k+q ) S(u-q + uj k+g )} 



(5.3) 



As in the scalar case, we are mainly interested in the infrared effects associated with 
the emission and absorption of soft photons in the intermediate state and only the finite 
temperature contribution. Therefore we will: i) neglect the contribution of the distribution 
function of the intermediate charged scalar, i.e. the term n k+q , ii) replace N q « in the 
expression and neglect the vacuum contribution (the one), thus obtaining 



p{k;u) 



e 2 kT 
' 2tt 2 
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with q c an upper momentum cutoff q c «T and 
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2(uj t k) 



The second contribution with support below the light cone is identified with the Lan- 
dau damping cut. Since the time dependent correlation functions involve the product of 
p(k; u>) cosout the frequency integral in the range — k < uj < arising from the Landau 
damping cut can be combined with the contribution from the positive frequency range by 



the symmetry of the integrand in eq. ( 3.20 ). After a straightforward calculation we finally 
find the finite temperature contribution to the spectral density near threshold to be given 
by 



Pi.R.(k;uj) 




nin 



UJ 



k 



uj + k 



+ 2Axu j [e(oj - uj k ) + &(k - uj)} 0(. 



(5.4) 



There are several noteworthy features of these spectral density near threshold, as com- 
pared to the simpler case of the scalar theory studied in section III (see eq. (fj.28|). In this 



case the spectral density is constant at threshold. The term uj 2 (multipling the logarithms in 
fl5.4p ) and the last term 8a; k arise from the region of the integral for which cos 2 9^1 which 
would give a vanishing integrand were it not for the fact that there is a linear divergence as 
uj — > k. These contributions arise from the emission and absorption of photons which are 
almost collinear with the incoming charged scalar. 
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We will now study the case of a massive scalar in two important limits: i) k/u k — v « 1 
{m ^ 0), ii) k/uf. = v m 1. In both cases the leading contribution is easily recognized to arise 
from the production cut oj > oj k . For m 7^ there are no infrared divergences associated 
with the Landau damping cut. 



1. v « 1 



Since the spectral density is slowly varying near threshold we can find the asymptotic be- 
havior in time of the coefficient of A k e tuJkt in (|3.20| ) at large times following the same method 
as in previous sections and appendix B. We find the following results: i) the imaginary part 
of the coefficient of A k e iuJkt is given by 
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(5.5) 



(5.6) 



the function < f(v) < 2 for < v < 1. 

ii) The real part of the coefficient of A k e tWkt in the asymptotically large time limit is 
given by (see appendix B) 
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+ O (e 2 T 
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(5.8) 



where we used the expression (|3.10|) and the fact that for the production cut, the spectral 
density only has support for positive frequencies. We then find the same phenomenon as in 
the previous case of the scalar particles in that the damping rate is one half of the expected 
value. The reason again is that the full spectral density for the charged scalar is very similar 
to that featured in figure ([!]) with a prominent peak near threshold that is almost half of a 
Breit-Wigner peak. However the logarithmic phase clearly exhibits the fact that cannot be 
interpreted clS 8b quasiparticle resonance. 
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The contributions from the coefficient of A* k e luJkt and the linear secular term from (/)% 
that contributes to the imaginary part are both subleading. We now renormalize the ampli- 
tude as in eqs. ( |3.23|) -( ^24]) with the choice 
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e 2 T 

— f(v) rln/ir 



TTUJk 



■ log(/i r e 1 



The solution of the dynamical renormalization group equation now leads to the asymptotic 
behavior of the expectation value of the charged fields 



(j) k (t) = A k (t ) e iVkit ' to) e - rk{t - to) 
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(5.9) 



Where we integrated the dynamical renormalization group equation with initial condition 
at to which is taken as some arbitrary renormalization point replacing the infrared cutoff /Z. 
The renormalization group invariance of 4>k(t) is now explicit, a change of the arbitrary scale 
to is compensated by a change in the amplitude Ak(to)- 



2. v 1 

The limit v — > 1 must be studied carefully. As uj — > Uk ~ k the term (oj 2 — k 2 ) In \u — k\ 
cannot be taken outside of the integral. However upon the change of variable uj — k = z/t 
in the integral in the same manner as that leading to ( |5.7| ) leads to a term that is of the 
form hx[fd]/t for this contribution, which then becomes subleading compared with the term 
in pi.R.(k; uj) that does not vanish as v — > 1. The asymptotic large time behavior is therefore 
obtained from the previous section with v ^ 1 by simply setting v = 1. The contribution 
to the final result in this limit arises solely from the emission and absorption of collinear 
photons. 

It is illuminating to try and understand this result in the hard limit v — > 1. The delta 
functions in ( |5.3| ) in the limit k ^> m become 8 (id — k — q =)= qcos8) therefore as uj — » k 
the whole contribution arises from photons that are emitted or absorbed collinearly 9 = 0, ir 
with the moving (hard) scalar. However, as pointed out originally by Pisarski |2U[| , the con- 



tribution from collinear photons does not survive screening effects arising from higher order 
contributions to the photon propagator. In particular, dynamical screening as a consequence 
of Landau damping of the intermediate photons cuts off the contribution of collinear pho- 
tons, and lead to a greater contribution of photons emitted or absorbed at right angles with 
respect to the moving charged particle. 
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The analysis of this section is illustrative of the power of the dynamical renormalization 
group to obtain the asymptotic long time behavior. For the case of QED, QCD or SQED, the 
analysis presented in this section in terms of the bare propagators for the scalars and photons 
has very limited validity. For soft external momentum, the infrared region of the internal loop 
requires HTL resummation of the internal lines and vertices [ PU| , P2"P§| . The main purpose 
of our analysis in this section, however, was to illustrate how the dynamical renormalization 
group is capable of revealing novel forms of relaxation with logarithmic corrections, power 
laws, anomalous dimensions etc. These alternative forms of relaxation cannot be found 
by attempting to describe exponential relaxation and computing an imaginary part of the 
self-energy on shell. 

We now include screening corrections via HTL resummation of internal lines. We will 
consider the case of hard external momentum for which only the internal photon line must 
be HTL resummed (the scalar is massive and hard). 



B. Hard thermal loop-resummed photon propagators 



We now focus on the case of hard external momentum of the charged scalar. In this case 
only the internal photon line receives HTL corrections, since the scalar in the loop is massive 
and hard, the vertex does not require resummation because one of the momenta into the 
vertex is hard (the scalar) [^,^,^]. Hence this situation is simpler than the case of soft 
external momentum that will be studied elsewhere. 

In order to include the leading order screening effects in the photon propagator, we 



must use the hard thermal loop resummed propagators [20]. The generalization of the HTL 
resummation program in the Matsubara formulation of finite temperature field theory to the 
real time formulation is described in detail in appendix A, we collect here only the main 
ingredients. 

The photon propagators can be written as in (|5.1| ) -(|5.2| ) but now with the resummed 
Wightmann functions (see appendix) 

Gt(t-1?) = Jdq p T (q ,q) [1 + N(q )} e*»W 
Qf{t ~ t') = J dq p T (q , q) N(q ) e'^^ 
where in the hard thermal loop limit the spectral density for transverse photons is given by 
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where u p (q) is the plasmon pole and Z(q) its (momentum dependent) residue, which will not 
be relevant for the following discussion. Inserting these propagators in the expression for the 
self energy, keeping only the term N(qo) ~ T/qo in the resulting expressions, and focusing 
on the hard scalar limit k m T ^> m we find that the self energy can be written in the form 
of a dispersion relation just as in eq. Q3.3p but with 

^ Tk f „2^„ f 1 JV n V 2\ /"j„ Pr(go,g) 

The infrared region corresponds to small q, and we find that for q « eT the integrand 
/?t(<Zq> q) I '% is strongly peaked at q = 0. Fig. 2 shows f>T(qo,q)q 2 /qo vs - Qo m units of 
eT j y/12. We find that for q « eT the photon spectral density is well approximated by |32| 

p T (g ,g) i r 



p(k; w) = / g 2 dg f rfX (1 - X 2 ) / rfg 5(u - k - q - qX) 

H z J J -l J On 



g Trg 2 T 2 + q 2 

r = g 3 



7ie 2 T 2 

When this Lorentzian distribution is integrated with smooth functions, it can be expanded 
in the width obtaining 

;pq^ "<(•)- £*"(»> + ■•■ 

and the infrared behavior is dominated by the S(qo). Finally we find the infrared behavior 
of the spectral density of the self-energy to be given by 

e 2 kT ri* dq ( ' (u - k) 2 \ e 2 kT \u - k\ 

p(k; u) w — / — 1 - ~ — 5~ ln 5 " 10 ) 

7r z J|w-fc| g \ q z J 7i z fi 

where q* ~ /x < eT is an arbitrary upper momentum cutoff, which physically is of the 
order of the plasma frequency. We thus see that dynamical screening originating in Landau 
damping for the photon propagator suppresses the contribution of collinear photons and as 
uo — > k the contribution to the self energy arises primarily from photons emitted or absorbed 
at right angles pO| . This is the same situation as in QED |32fl . 

The final spectral density for the self energy, eq. ( |5.10| ) is therefore of the same form as 
the finite temperature self energy of the simple scalar theory studied in the earlier sections, 
see eq. ( |3.28| ), and thus justifies our excursion into that simpler theory. 

We can now follow the same steps to study the secular terms in the real time perturbative 
expansion, which lead to eq. ( |3.29j ). 

The long time asymptotic behavior is obtained by inserting the spectral density eq. 
( |5.10| ) in the expressions (|3.20|) - (|3.21|) with u k = k. The asymptotic dependence on time is 
extracted by changing variables too; — k = z/t and upon setting A k = A* k , we find: i) the 
imaginary contribution to (fp-' a ) + (f)( l,b > is given by an infrared finite and linear in time secular 
term which is interpreted as a renormalization of the frequency, ii) the real contribution to 
(j)(X> a ) -)_ is given by 

e 2 

Re<50 fc (t) = — aTt ln/it ; a — — 

47T 
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and the first order correction is thus found to be given by 

^^{t) — A k e Wkt [i 5uo k t — a T t ln/it] + c.c + regular perturbative terms 
Introducing the renormalization of the amplitude as in equations (|3.23|) - (|3.24|) and choosing 

\zr (t) = aT t hx/jLT ; \zr (t) = —5u k r 
we obtain the renormalization group equation 

- [i6u k - a T (In /it + 1)] A k (r) = 

with solution 

A k (r)=A k (t ) e ^(-*o) e -" T ^r 

where we have chosen to = fi^ 1 as an initial condition for the integration. Now choosing the 
arbitrary renormalization scale r to coincide with the time t we finally arrive at one of the 
main results of this article which is the asymptotic behavior in time of the expectation value 
of the charged scalar field 

(f) k (t) = A k (t ) e ^(*-*o) e -«"M*/*o] + c.c. (5.11) 

thus displaying the renormalization group invariance of the solution: a change of the arbitrary 
time scale to (inverse of the infrared cutoff) is compensated by a change in the amplitude. 
The quantity u k is the renormalized mass including the (infrared finite) HTL corrections. 

A similar behavior for the asymptotic dynamics of the fermion field in QED has been 
obtained via the Bloch-Nordsieck resummation in reference ||32|| . The power of the dynamical 
renormalization group has now become explicit in that it transcends any approximation and 
implements a resummation of the logarithmic infrared divergences that in this case includes 
the resummation of the hard thermal loops. 



VI. CONCLUSIONS, COMMENTS AND MORE QUESTIONS: 

In this article we have introduced a novel method of dynamical renormalization group 
resummation to study relaxation in real time. The first step of the program is to relate the 
retarded Green's functions that contain the dynamical information on the time evolution 
away from equilibrium in linear response to an initial value problem for the expectation 
value of the fields. This initial value problem in real time allows to implement the method 
of dynamical renormalization improving the perturbative solution by a resummation. This 
resummation is the real time counterpart of the resummation via the renormalization group 
in Euclidean field theory. We first apply our methods to a scalar theory of one massive and 
one massless field. The emission and absorption of the massless field introduced infrared 
divergences akin to those found in gauge theories. We compare in such a model the late 
behaviour of the massive field amplitude at zero temperature |eq. (|3.27|) ] to results based on 
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Bloch-Nordsieck resummations [eg. ( |3.16| )1 as well as the renormalization group applied to 
Euclidean Green's functions [eg. (|3.18|) 1. 

Furthermore, for large temperature we compute the the late behaviour of the massive field 
amplitude both for hard and soft external momenta. For hard modes, the field amplitude 
relaxes as e - *- 9 T / 47r )' M*/*o] where g stands for the coupling constant. 

In real time the infrared divergences are manifest as secular terms in the perturbative 
solution of the evolution eguation for the expectation value of the fields. The dynamical 
renormalization group implements a resummation of this secular terms that leads to an 
asymptotically convergent solution and clearly describes relaxation in real time. 

After applying our method to the scalar model, we focused our attention to implement- 
ing the dynamical renormalization group to resum the infrared divergences associated with 
massless transverse photons in scalar QED at finite temperature. The infrared divergences 
in this theory are similar to those found in QED and in lowest order in QCD. 

We have included the resummation of the hard thermal loops and Landau damping in 
the internal transverse photon propagators, and implemented a dynamical renormalization 
group resummation. The renormalization group improvement leads to an anomalous loga- 
rithmic relaxation for hard modes as a conseguence of infrared divergences associated with 
the emission and absorption of photons at right angles. These anomalous logarithmic relax- 
ation are similar to the scalar field behaviour and consistent with those found in QED via 
the Bloch-Nordsieck resummation f32fl . 

In all cases investigated (both in the scalar model and in QED) the field behaviour 
prevents an interpretation of the relaxation of charged excitations in the medium in the 
form of a simple exponential with a damping rate determined by the imaginary part of the 
self-energy on-shell. 

The advantage of the dynamical renormalization group is that its implementation is 
rather simple and trascends any approximations of the Bloch-Nordsieck type, it can be 
consistently improved by considering higher orders in the hierarchy of eguations obtained in 
perturbation theory. 

Furthermore, the real time dynamics obtained via this resummation program leads to a 
clear interpretation of the relaxational processes and time scales without any assumptions 
on the validity of the guasiparticle picture of collective excitations. The analysis of secular 
terms in lowest order provides a simple criterion for deciding if the collective excitations can 
be described as narrow resonances with a width determined by the imaginary part of the 
self-energy on shell: linear secular terms lead to such a guasiparticle description, non-linear 
secular terms in lowest order signal anomalous, non-exponential relaxation. 

The dynamical renormalization is a different resummation scheme than the HTL re- 
summation, and the latter can be consistently included in the former as was shown in this 
article. 

We are currently implementing this method to study relaxation of soft and hard fermion 
and gauge fields in QED and QCD directly in real-time, thus bypassing the conceptual 
limitations of the guasiparticle picture in the sense of exponential relaxation and a damping 
rate determined by the imaginary part of the self-energy on-shell. We expect to report results 
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in the near future. 

Comments and further questions: As we have seen in the example worked out 
in detail in section III, at finite temperature the infrared divergences are akin to those of 
a superrenormalizable theory of critical phenomena, in the sense that they are no longer 
logarithmic because the temperature introduces a new scale. There are very few methods 
for renormalization of infrared divergences in superrenormalizable theories near the critical 
point, one of the most popular being the e expansion where e = d — 4 is the departure from 
the upper critical dimension. Whereas the validity of the epsilon expansion appended by 
Pade resummation has been confirmed in Ising-like models via either strong coupling lattice 
expansion or Montecarlo simulations, the validity in a general case is at best questionable 
for e = 1. The e expansion resums some subset of the Feynman diagrams and only some 
part of them E71 the leading logarithms. 

The dynamical renormalization group therefore provides an alternative to study the in- 
frared divergences directly in real time by resumming secular terms in the perturbative 
solution of the equation of motion. In the case of infrared divergences the secular terms 
reflect these in the form of logarithmic dependence on time. However, the uselfullness of 
the dynamical renormalization group is not restricted to this logarithmic divergences, as 
explicitly shown in section III, in the usual case of narrow resonances, the secular terms are 
linear (in lowest order) and their resummation through the renormalization group equation 
leads to the usual quasiparticle real-time evolution. 

There is a translation between the resummation implied by the usual (Euclidean) renor- 
malization group and that by the dynamical version: the Euclidean version sums the leading 
logarithms |36| , f4"T|| the dynamical version sums the leading secular terms [B7 



It has been shown in ref. [0 using the formal theory of envelopes that the dynamical 



renormalization group resummation of secular terms provides an uniform approximation 
to the exact solution for systems of ordinary differential equations. This is true to any 



given order of perturbation for arbitrary ordinary differential equations |37|- It will be very 



interesting to extend such a proof to the evolution equations considered in the present paper. 

Thus, for the moment, the situation with the dynamical RG is similar to that of the e 
expansion in critical phenomena for e = 1: it provides a resummation scheme for the infrared 
behavior in a consistent manner and it agrees with known results in cases where it can be 
compared. Furthermore in the case in which the renormalization group and Bloch Nordsieck 
lead to non-trivial exponentiation of infrared divergences, the dynamical RG reproduces the 
results in real time. Thus we believe that the cases analyzed in detail in this article and those 
analyzed in the literature provide very strong evidence for the validity of this approach. The 
promise of the dynamical RG as a powerful method to study transport phenomena warrants 
a deeper study on the renormalization aspects of the evolution equations in real time and a 
more formal proof of the applicability of the dynamical RG in these problems. This avenue 
of study is currently in progress. 
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APPENDIX A: EXACT RETARDED PROPAGATORS 

In this appendix we gather and generalize some results of the HTL resummation pro- 
gramme in Matsubara finite temperature field theory |1(J, to real time. 



Consider the exact equilibrium Wightman and retarded Green's functions for a real scalar 
field 



- iG>(t - = = -ij du G>{u) e-^'> 
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J ATX 
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Inserting a complete set of eigenstates of the full interacting Hamiltonian, we obtain the 
spectral representations for the Fourier transforms given by 

- iG>(u) = ±Y, e ~ PEm IH$fc(0)H| 2 ^ - (E n - E m )) 
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where the last equality, the KMS condition is obtained by relabelling m —>■ n in the sum and 
Z = Em e~ f3Em is the equilibrium partition function. Inserting a representation of the theta 
function we finally obtain 



GR,k(qo) = - / du- 
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p(u, k) = -i \G>{u) - G<{uj)} = -Im[G R>k (u)} 
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Using the KMS condition ( |A1| ) we can finally write the exact non-equilibrium Wightmann 
functions in terms of the exact spectral density as follows 

= fdq p(qo,k) [1 + N(q )} 
(*-S(f )*£(*)> = I dQo p(qo,k) N(q ) e**W 
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Using the KMS condition, and relabelling the sum indices in the spectral representation 
we find that p(qo, k) = —p(—qo, A;). The same steps lead to an equivalent expression for the 
transverse gauge fields, whose Wightmann and Green's functions are proportional to the 
transverse projection operator. 

The advantage of the representation ( |A2[ ) is that once we compute the spectral function 
p(qo, k) in some approximation, we can insert the Wightmann functions in the internal 
loops thus providing a resummation of the perturbative series. The results from the first 
section allow to obtain the retarded propagator Gr^Qo) from the solution of the initial value 
problem with an external source through the relation to linear response as detailed in the 
first section. For example by studying the equation of evolution for the expectation value 
of the transverse photon fields in the HTL approximation as was done in [45] we can obtain 
the spectral representation of the transverse fields in the HTL approximation and the results 
of this appendix allow us to implement a resummation of screened photon propagators into 
the real time description. 



APPENDIX B: ASYMPTOTIC BEHAVIOUR OF SPECTRAL INTEGRALS 

We summarize in this appendix the late time behaviour of integrals over the density of 
states used in Sections III and IV. 

In the formulas below p(y) stands for a smooth function for < y < oo. p(0) as well as 
p'(0) are finite. For large y ,p(y) decreases as a power such that the integrals over y converge 
at infinity. These properties are fullfilled in all cases where these formulae were used in the 
paper. 

/ 4 (1 - cos yt) p(y) =°° - t p(0) + p'(0) [ln(/x t) + j) 
Jo y 2 

+ l Q X ^\p(y)-p(0)-yp'(0)e(p-y)} + o(^j , 

r - C 1 - cos ^) p(y) * =°° p(o) Ka*«) + 7] + r - \p(v) - P(0) 0{n -y)] + o (I) , 

Jo y Jo y \t / 

r^(t- S ^y±) p( y ) ^ t p(0) pn^i i)+7-l]+if- \p(y) - P(0) 0(p-y)] + O (±) , 
Jo y \ y J Jo y \tj 

where 7 = 0.5772157. . . is Euler's constant. 

Notice that the formulas are independent of the scale p as one can easily see since 
the derivative with respect to fi of the r. h. s. identically vanishes. The scale p has been 
introduced just to have a dimensionless argument in the logs. 

We have also used similar integrals when the resonance was away from the threshold. We 
have for such case, 

r % (1 - cosyt) p(y) 1 =°° vr t p(0) + V r % \p(y) - P(0)] + O (\) , 
J -A y z J -A y z \tj 
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-a y 



y 



where A is a fixed positive number. 

In addition, we needed in sec. Ill and IV integrals where the spectral density has a 
logarithmic singularity at a finite point. 



dy 



^ (1 - cos yt) p(y) In M * =°° tt t p(0) 1 - 7 - ln(2 * T) 



2T 



A 1/ 



2f 



+ J>£ *[*,)-,*»] lnM + (I) , 



1/ 

' - l°° — p(y) lnM-Trp'(O) [ln(2tf)+ 7 
Ay 2T l 



+ 



APPENDIX C: A SIMPLE EXAMPLE: THE DAMPED HARMONIC 

OSCILLATOR 

In this appendix we give a rather simple example of the dynamical renormalization group 
for pedagogical reasons and to illustrate the fundamental features within a simple setting. 
We consider the equation of motion of a damped harmonic oscillator: 



y + y = -ey,e«l 



(CI) 



and seek a solution in a perturbative expansion in e of the form y = y + ey\ + e 2 y2 + 
where the yi are solutions to the following hierarchy of equations: 

Do + Vo = , 
Vi + Vi = -2/o , 

$2 + yi = -iii , 



These equations can be solved iteratively by starting from the zero order solution 

y (t) = A e u + c.c , 
in terms of the retarded Green's function 

G iet (t - t') = sm(t - t') 6{t-t') , 
Up to second order in e, the solution is given by 



y(t) = A e 



it 



+ c.c. + non-secular 
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Note that this solution contains secular terms that grow in t, the terms denoted by non- 
secular remain finite at all times. We see that the perturbative expansion breaks down at 
a time scale « 1/e. The expression in the brackets can be interpreted as a change in the 
complex amplitude. The dynamical renormalization is achieved by introducing a time scale 
t at which the secular terms are absorbed in a renormalization of the complex amplitude. 
We write A = A(t) Z(t) with Z(f) = 1 + e Z\{r) + e 2 z 2 (r) and choose Zi(r) to cancel the 
secular terms at the scale r, this is similar to choosing the renormalization scale in the usual 
renormalization program. Up to 0(e 2 ) we find 

t r 2 t 
Zi(t) = — ; zo(r) = i — . 

IK J 2 , A J g g 

After renormalization the solution is given by 



y(t,r) = A(r) 



1 



2 v 



e 2 e 2 
r) + — (t - r) 2 + i — (t 



T 



+ c.c. + nonsecular . (C2) 



Since r is an arbitrary scale, the solution cannot depend on it, thus the statement 
dy(t, T)/dr — leads to the dynamical renormalization group equation to this order 



dr K J \ 2 8 . 







where we have expanded the r derivative of the amplitude in a power series expansion in e 
consistently to second order. Obviously the solution to the renormalization group equation 
is given by 

A(r) = A(0) e"5 T e 1 ^ 
setting t = t in QU2p we finally find 



y(t) = A(0) e 



fr e i(l-V)* 



+ C.C. 



which is obviously the correct solution to second order. Further simple and not-so-simple 



examples can be found in ref. [37 
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FIGURES 




FIG. 1. The spectral density for the a field S(k,uj;T) (eq. ( 3.33[) ) near threshold for 
p(k,io k ) = -0.005. 
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FIG. 2. p T (q , q)(q z /qo) for q = 0.1^= vs. q in units of ^= 



eT 
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